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Abstract 

Finite 3D Ising ferromagnets are studied in periodic magnetic fields both by 
computer simulations and mean-field theoretical approaches. The phenomenon of 
stochastic resonance is revealed. The characteristic peak obtained for the correlation 
function between the external oscillating magnetic field and magnetization versus 
the temperature of the system, is studied for various external fields and lattice sizes. 
Excellent agreement between simulation and theoretical results are obtained. 

PACS number(s): 05.40. +j; 75.10.-b; 02.70.Lq 



* Permanent address: Babe§-Bolyai University , Dept. of Physics, str. Kogalniceanu 1, RO-3^00 
Cluj-Napoca, Romania 



1 



1 Introduction 



It is well-known [1-4], that periodically modulated bistable systems in the pres- 
ence of noise exhibits the phenomenon of stochastic resonance (SR). The bases of 
this phenomenon is that the correlation a between the modulation signal and the 
response of the system presents an extremum for a given noise intensity. 

Evidences for SR were found in analog simulations with proper electronic circuits 
[3,5,6], laser systems operating in multistable conditions [7], electron paramagnetic 
resonance [8,9], in a free standing magnetoelastic ribbon [10] or in globally coupled 
two-state systems [11,12]. Theoretical aspects of the problem were reviewed in [4]. 

Considering a special and practically important case of coupled two-state sys- 
tems, recently we reported on the possibility of obtaining the SR in finite two dimen- 
sional Ising systems [13]. In the mentioned work, in contrast with all other earlier 
works we did not consider any external stochastic forces, only the thermal fluctua- 
tions in the system. The problem was studied by computer simulations considering 
a heat-bath dynamics. 

The present paper intends to complete the earlier one [13], studying the impor- 
tant three dimensional (3D) case. Presuming that in 3D the mean-field theories give 
reasonable results, we also consider a theoretical approach to the phenomenon. We 
discuss and compare our theoretical and computer simulation results. 

2 The Method 

For detecting SR we need a system in a double well potential, governed by a stochas- 
tic force. In the meantime the two minima of the double well potential must be mod- 
ulated periodically in antiphase. One can immediately recognize that ferromagnetic 
Ising systems in oscillating magnetic fields satisfy all these conditions: 

• at zero thermodynamic temperature the free-energy versus magnetization curve 
has the double well form. 
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• an external periodic magnetic field [B{t)], modulates the two minima in an- 



• the effect of a positive temperature can be viewed as a stochastic driving force, 

• the magnetization as a function of time {M{t) = J2iSf) can be considered as 
the response function of the system. 

Due to the fact that in the proposed system the noise intensity (thermal fluctuations) 
is temperature dependent, the characteristic maximum for SR must be detected at 
a resonance temperature (T^) in the plot of the correlation a {a —\< B{t) ■ M{t) >\) 
versus the temperature (T) . 

The Hamiltonian of our problem is 



where the sum is referring to all nearest neighbours, Sf — ±1, hb is the Bohr 
magneton and B{t) the external magnetic field. We consider B{t) in a harmonic 
form: 



We mention that recently 3D Ising systems in oscillating magnetic fields were 
already considered by computer simulations [14] . However in [14] the authors study 
the hysteretic response of the system and no evidences for SR are discussed. 

We will study the proposed system (1) from the viewpoint of SR both by com- 
puter simulations and by a theoretical approach. 

2.1 Computer Simulation 

To study the time evolution of the proposed system (1) first we considered a com- 
puter simulation with heath-bath dynamics. The time scale was chosen in a conve- 
nient way, setting the unit-time interval equal with the average characteristic time 
(r) necessary for the flip of a spin. We have taken this time interval r as con- 
stant, and thus independent of the temperature. Although this assumption is just a 



tiphase. 



H^-JY^StS| + ^JiBB{t)Y,St, 



(1) 




(2) 
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working hypothesis we expect to give useful quahtative results. The spin flips were 
reahzed with the probabilities of the Metropolis [15] algorithm, choosing the spins 
randomly at each moment. 

The simulations were performed on cubic lattices with W = N x N x N spins, 
considering the value of N up to 50. One simulation step was defined as W trials 
of changing spin orientations, and corresponds to a time interval r. The period (P) 
of the oscillating magnetic field will be also given in these r units. The amphtude 
A of the magnetic field is considered already multiplied by HB/k, and thus will 
have the dimensionality of the temperature {k is the Boltzmann constant). The 
temperature will be considered in arbitrary units, and the critical temperature of 
the infinite system {Tc 4.444 J/k) will be set to 100 units. Starting the system 
from a random configuration we considered 10000 simulation steps to approach the 
dynamic equilibrium. The correlation function between the driving field B{t) and 
the magnetic response M(t) 

1 " 

a ^\< B{t) ■ M{t) >\ = \ -Y,B{U)M{U) \ (3) 

^ i=i 

was studied after this during 10000 extra iterations. (The averaging in (3) is as a 
function of time, and ti — r-i). 

The correlation [a) was studied function of: 

• the temperature (T) , 

• the lattice size (N), 

• the amplitude of the magnetic field (A) , and 

• the period of the oscillating magnetic field (P). 

2.2 Theoretical approach 

We propose now a mean-field like theory to describe the time evolution of the system. 

Let us consider that P{M, t) dM is the probability for having the magnetization 
of the system between the values M and M + dM at a time moment t. We know 
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from the classical fluctuation theories that for a system in thermal equilibrium with 
a heat-bath the probability P{x) of getting into a state x is proportional with the 

P{x) oc e-^-^(^) (4) 

factor, where [3=1/ (kT), and F{x) is the free-energy in state x. Presuming that for 
a very short time interval (dt) and fixed Mq the changes in P{Mo, t) are influenced 
mainly by the neighbouring values P{Mo + dMo,t) and P{Mo — dMo,t), the time 
evolution ol P{M,t) can be approximated as 

3P(M f) (p-l3F{M+dM,t) , -f3F{M-dM,t)\ 

-^^^ - ^ ^-P{M^ t) + (5) 

-^F(M,t) 

+ ^ [P{M - dM, t) + P{M + dM, t)], 

where: 

S=Y1 e-^^(^'*). (6) 

{M} 

Up to second order terms one can consider: 

P^M + ,M.t) = PiM,t) + 'J^,M + l?^,M^ (7) 
P(M - ,M, t) ^ nu, t) - '-^,M . 

In this manner equation (5) will be written as 

dP{M,t) d'PjM^t) d'FjM^t) aF(M,t) 

-^r~ - dM^ + dM^ - ^ ^-^Ar~^ (s) 

with: 

-/3F(M,t) 

C = ^— dW. (9) 

A few remarks on the obtained approximative "master" equation (8) must be 
made: 

• The equation is reasonable from the point of view that for dF{M,t)/dt = 0, 
P(M, t) = Const, exp [— /3F(M, t)] is a stationary solution. 



• The equation will not necessary keep the normalization of P{M, t). For assur- 
ing that P{M,t) will have the meaning of probability, one must apply equation 
(8) together with the 
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w 

P{M,t)dM = l (10) 

-W 

normalization condition at each time step. {W is the number of spins in the 
system.) 

• The value of C determines the speed of changing P(M, t) at each moment. 
In reality C is temperature dependent, but as a first approximation one can 
consider it constant. Fixing it at a reasonable value we assume that will not 
influence a quahtative discussion on the phenomenon of SR, only the absolute 
values of the a correlation. 

The F{M,t) free energy of our system of W coupled = ±1 spins in a B{t) 
magnetic field can be computed using a mean-field approximation. We assume that 

F{M,t) = U{M) -T.S{M), (11) 

with U the internal energy, T the temperature and S the entropy. For a cubic 
system, in the mean-field approximation we consider 

U{M) = -3J^ + I^bM ■ B{t), (12) 

and: 

( W+M \ 
2 ). (13) 
w I 

The first and second order derivatives of F{M,t) can be approximated numeri- 
cally. Equation (8) with the condition (10) can be now solved also numerically for 
a given C value. Starting from a P{M, 0) initial distribution in principle we should 
be able to compute P(M, t). The a correlation can be calculated as 

a^- dt M- B{t) ■ P{M, t) dM, (14) 
r Jo J-w 

which is also computable. 



Choosing the same units for J and B as in computer simulations we did solved 
numerically equation (8) with the normalization condition (10). B{t) was taken in 
the form given by (2) and we fixed J — 20. The "critical temperature" of the system 
without magnetic field became in this way T^. 122. We also fixed W = 4Q and 
considered different temperatures (T) around T^, several periods (P) and several 
amphtudes (A) for the oscillating magnetic field. The constant C dt was taken as 
0.01. We also checked that the overall picture of a versus T is not significantly 
infiuenced by lowering the value oi C dt. 

Unfortunately the time scale of this theoretical approximation can not be obvi- 
ously related to the time scale used in computer simulations. Comparision between 
theoretical and computer simulation results must be viewed under this assumption. 

3 Results 

Our simulation results and theoretical approximations are summarized in Figs. 1-6. 

In Fig.l we present a characteristic computer simulation result for the shape 
of the (7 versus T curve. One will observe that in accordance with the predicted 
phenomenon of SR, at a given temperature a presents a maximum. The tail of 
this resonance peak is nicely described by a power law (bottom picture) . Analysing 
the scaling exponents for different simulations (different P, A and values) wc con- 
cluded that it lies in the [(-1.7)-(-2.4)] interval. A characteristic theoretical picture 
for the same a versus T curve is presented in Fig. 2. Again the resonance peak is 
observed and nice concordance with simulation results are suggested. Moreover the 
tail of the obtained resonance peaks are also described by power laws with scahng 
exponents in the [(-1.94)-(-2.05)] interval. 

Figs. 3 and 4 presents the shape of the resonance peak for several values of the 
modulation amplitude A. Both theory (Fig. 4) and simulations (Fig. 3) predicts 
that the resonance temperature is almost independent of the modulation amplitude, 
exhibiting only very slight variation as a function of this (i.e. for higher amplitudes 
Tr is shifted in the direction of smaller values). In contrast with this, the height of 
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the peak depends sensitively on the values of the modulation amplitude. 

As we concluded also in [13], for small lattices the Tj. resonance temperature is 
strongly dependent on the lattice size (N), and in the limit of relatively big lattices 
[N K, 20) is tending to a constant limiting value. In this sense for three different 
modulation periods P, the simulation results are presented in Fig. 5. We did 
not studied this dependence theoretically because our numerical calculations were 
technically hmited by the values of W (i.e. only small lattices with 1^ < 50 were 
computable). Prom Fig. 5 we also learn that the temperature % is dependent on 
the modulation period P. 

We studied the dependence of T^. versus P both by simulations and theoretically. 
The results in this sense are plotted in Fig. 6. Both simulations and theory predicts 
that for high periods the % resonance temperature is tending to the Tc "critical 
temperature" of the system. This dependence can be described by a 

Tr-T^^ K^e-^^^ (15) 

exponential law. Both for simulations and theoretical approximations the step-like 
form of the results are due to the fact that in detecting Tj. the temperature was 
varied with steps of 5 units. The differences between the obtained K2 values in 
simulations and theory are due to the already mentioned problem that the time 
scales are not related each to the other. 

4 Conclusions 

The first conclusion would be that both our computer simulations and theoretical 
data suggest that the phenomenon of SR should be detected when one studies finite 
ferromagnetic systems in oscillating magnetic fields. The characteristic peak of SR 
is obtained by studying the a —\< B{t) ■ M{t) >\ correlation as a function of 
temperature. For a given resonance temperature T^, this correlation a exhibits a 
maximum. 
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Fixing the frequency, for small lattices {W < 4000) the Tr resonance tempera- 
ture is dependent on the lattice sizes and is tending to a limiting value for relatively 
large (W > 4000) lattices. The resonance temperature proved to be dependent also 
on the period of the magnetic field, and in the limit of large periods is converging 
exponentially to the critical temperature of the system. Because in real experimen- 
tal conditions we are in the very high period limit (the time unit in simulations is set 
by T, the characteristic time for the flip of a spin), we expect Tr to be detected at Tc. 
We also concluded that the Tr resonance temperature is not significantly influenced 
by the amplitude A of the oscillating magnetic field, the value A determining mainly 
the height of the resonance peak. 

Our theoretical results proved to be in excellent concordance with the computer 
simulations, justifying thus the considered approximations. 

We consider that from the viewpoint of magnetism an experimental study of the 
problem would also be of interest. 
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Figure Captions 



Fig. 1. Characteristic peak for tlie SR phenomenon obtained by computer 
simulations. The bottom picture illustrates the |< B{t) ■ M{t) >|= 11.181 • T"^ '''^ 
power law behaviour of the tail (m = M/W, N ^ 20, A ^ AO, P ^ 100 and 
T, = 100). 

Fig. 2. Characteristic shape of the SR peak given by our theoretical approxi- 
mation. The bottom picture illustrates the |< B(t) ■ m{t) >\= 2.088 • T~^'^^ power 
law behaviour of the tail (m = M/W, W ^ 40, J ^ 20, % = 122, P = 800 and 
A^l). 

Fig. 3. Computer simulation results for the shape of the SR peak considering 

several values of the modulation amplitude A {m = M/W , N = 20, = 100 and 
P = 50). 

Fig. 4. Theoretical results for the shape of the resonance peak considering 
several values of the modulation amplitude A {m — M/W, 1^ = 40, J = 20, 
rc = 122andP = 800). 

Fig. 5. Computer simulation results for the dependence of the resonance 
temperature (T^) versus the lattice size {N). Results for three different modulation 
periods (P) are presented (Tc = 100 and A — 10). 

Fig. 6. Dependence of the resonance temperature [\n{Tr—Tc)] versus the 
modulation period (P). The best fit line indicates the — Tc = 4.061 • e~° °^^ and 
T^ — Tc = 4.312 -6"° °°^^ dependence for simulation and theoretical data respectively 
( for simulation: A = 10, N = 20; and for theory: A = 1, W ^ 40). 
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